Thawing models in the presence of a Generalized Chaplygin Gas 
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In this paper we consider a cosmological model whose main components are a scalar field and 
a Generalized Chaplygin Gas. We obtain an exact solution for a flat arbitrary potential. This 
solution have the right dust limit when the Chaplygin parameter A — > 0. We use the dynamical 
systems approach in order to describe the cosmological evolution of the mixture for an exponential 
self-interacting scalar field potential. We study the scalar field with an arbitrary self-interacting 
potential using the "Method of /-devisers". Our results are illustrated for the special case of a 
cosh-like potential. We find that usual scalar field dominated and the scaling solutions cannot be 
late time attractors in the presence of the Chaplygin gas (with a > 0). We recover the standard 
results at the dust limit (A—tO). In particular, for the exponential potential, the late-time attractor 
is a pure generalized Chaplygin solution mimicking an effective cosmological constant. In the case 
of arbitrary potentials, the late-time attractors are de Sitter solutions in the form of a cosmological 
constant, a pure generalized Chaplygin solution or a continuum of solutions, when the scalar field 
and the Chaplygin gas densities are of the same orders of magnitude. The different situations depend 
on the parameter choices. 

PACS numbers: 98.80.Cq 



Introduction 



It is of common knowledge that the expansion of 
the universe is currently passing through an accelerated 
phase p], • All the observational data from these for- 
mer references until the current measurements of red- 
shift and luminosity-distance relations of Type la Su- 
pernovae (SNe)j;3j are in agreement with this accelerated 
era. These observations are indicating the presence of 
a vacuum energy, the well known cosmological constant 
Another possible description is the existence of 
the scalar field that is evolving in a universe described 
by a Friedmann-Robertson- Walker (FRW) geometry, the 
quintessence model |7J ■ The current acceleration of the 
universe is one of the important problems of the mod- 
ern cosmology. This problem appears in the Einstein's 
standard general relativity, and one of the porposals to 
solve it, within this framework, is to consider one ex- 
otic component in the matter content of the universe Q , 
the so called dark energy component (for a review about 
this see Ref. d). We would like t0 focus 

our attention 

in models where dark energy is described by a homoge- 
neous scalar field (</>(£)) In order to describe the 
dynamics of dark energy (quintessence) the equation of 
state parameter (EoS), uj^ = p^/ ' plays a crucial role. 
Its current value is close to — 1. Therefore the cosmic 



evolution of dark energy through the scalar field is de- 
scribed by the barotropic parameter ujq. In Ref. fTTI j . 
were distinguished between two categories, the freezing 
(dujc/j/dcj) < 0) and thawing (dw^/dcj) > 0) models. These 
quintessence models are characterized for a, scalar field 
potential that asymptotically goes to zero pj|. There are 
several references about this kind of cosmological mod- 
els, for a summary see the Refs. p~3|. On the other hand, 
Chaplygin gas models have been widely investigated in 
the literature [TH-H?! . The Generalized Chaplygin Gas 
(GCG) have been investigated from the dynamical sys- 
tems viewpoint for example in |28-(32j . 
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In this article we would like to extend the analysis in 
[l2| by considering a more general matter component, 
that is a GCG and a scalar field characterized by its 
self-interacting scalar potential. We obtain an exact so- 
lution for a flat arbitrary potential, that have the right 
dust limit when the Chaplygin parameter A — > pl| . In 
order to motivate the analysis for a general (arbitrary) 
potential, we first consider the simple case of an exponen- 
tial potential V{4>) = Voe~ A( ^ [33l - t52 |. Then, we study an 
arbitrary self-interacting scalar field potential. In this 
case we use the "Method of /-devisers" presented in [53[ . 
This method allows to perform the phase space analysis 
without specifying the potentials ad initium, and then 
one just substitute the desired forms, instead of repeat- 
ing the whole procedure for every distinct potential. The 
method is a refinement of a method that has been ap- 
plied to isotropic (FRW) scenarios [H, lo4l - [58| , and that 
has been generalized to several cosmological contexts [59i — 

m. 

This article is organized as follow, in section U we 
present the cosmological model under consideration. Sec- 
tion [XT] is devoted to the study of the exponential poten- 



2 



tial an the corresponding dynamical system. In section 
IIIII we study the phase space of the cosmological model 
for general potential V(4>). Finally, we conclude in the 
section [TVl 



I. THE COSMOLOGICAL MODEL 

In this section we would like to describe a cosmologi- 
cal setting compose by a minimally-coupled scalar field 
that describes the dark energy and a Chaplygin gas that 
behaves as dark matter in the appropriate limit. 

The cosmological equations are given by 



„2 Pch + P$ 
H 3— 

p' ch + 3H(p ch + P ch ) 
P4, + 3iJ(p + P ) 



= 0, 
= 0, 

= 0, 



(1) 

(2) 
(3) 



where we work in units in which 8irG = 1, H is the Hub- 
ble constant, p c h and p^ are the Chaplygin gas and scalar 
field densities, respectively. For simplicity, the radiation 
component is neglected. P c u and P^ represent the Chap- 
lygin gas and scalar field pressures, respectively. For the 
scalar field we have: 



P/> 



On the other hand, the EoS for the GCG is 



P. 



ch 



A 

Pch 



(4) 
(5) 

(6) 



where A is a positive constant and a is a constant with 
an upper bound, a < 1. In particular, when a = 1 corre- 
sponds to the original Chaplygin gas. In the framework of 
FRW cosmology, this EoS leads, after inserted into the 
relativistic energy conservation equation, to an evolution 
of the energy density as 



and BOOMERANG (by using the third peak localiza- 
tion). Recently, the values B s = 0.7751°;°^^° °^ and 



+0.000970+0.00268 
-0.00126-0.00126 



were obtained from Markov 



a = 0.00126 

Chain Monte Carlo method |65( . For the numerical sim 



ulations implemented in the present paper we select the 
value a = 0.5. This value seems to be large in comparison 
with the above observational values, however, if we con- 
sider large values for a, and include the effect of shear 
and rotation, then when studying the evolution of the 
perturbations in GCG universes, it is found that that 
the joint effect of shear and rotation is that of slowing 
down the collapse with respect to the simple spherical 
collapse model. The described effect allows to solve the 
instability problems of the so-called unified dark matter 
models at the linear perturbation level (66[. 

The evolution of the energy density p c h shows the be- 
haviors of GCG at different times. At early times, the 
energy density behaves as matter while at late times it 
behaves like a cosmological constant. Then, this GCG in 
principle describes both dark matter and dark energy in 
a single matter component. 

Now we can define the barotropic index uj c h'. 



A 



A 



(8) 



At late times, the universe is dominated by the scalar 
field and the GCG, neglecting the radiation component. 
From the Friedmann equation (JI]) for a flat universe filled 
with a scalar field and GCG and defining the density 
parameter for each component as fix — 3^2 we found 
that our components satisfy il^ + Q c h = 1. Equations 
@ and ([3]) can be rewritten in terms of the auxiliary 
variables x, y and A, defined by 



V6' 



y = 



3H 2 ' 
1 dV 

V ~d~6 : 



(9) 



Pch 



A 



B 



7 3(a+l) 



PchO 



(1 - B 3 ) 



a 



3(0+1) 



Here, a is the scale factor and B is a positive integra- 
tion constant. In this way, the GCG is characterized 
by two parameters, B s = Aj p 1 ^ and a. Here, p c ho 
is the current value of p c h, considering that a = 1 at 
the present. These parameter has been confronted with 
the observational data, see Refs. [63l. |64| . In particu- 



0.73 



+0.06 
-0.06 



and 



-0.09 



+0.15 



-0.12 

the values 



lar, the values of B, 

were obtained in Ref.[64j. Also, in Ref.fl^ 
0-81 < B s < 0.85 and 0.2 < a < 0.6 were found from the 
observational data arising from different colaborations, 
such that Archeops (by using the first peak localization) 



where the prime denotes derivatives with respect to t = 
In a. Considering that the contribution for the kinetic 
and potential energy are given by x 2 and y 2 respectively, 
the density parameter of scalar field is given by 



3H 2 

therefore the equation of state is 

7 = 1 + LU = 



2 1 2 

x +y , 



2x 2 



(10) 



(11) 



x 2 + y 2 

Thus, inserting the auxiliary (JSJ into the equations of 
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motions (H),© and ([3]) we arrive to the following system 



-3x 



-sy 2 + ^{2x 2 + (1 



x 2 -y 2 )(l + uj ch )}, 
(12) 



^{2x 2 + {l-x 2 -y 2 ){l + L0 ch )} 1 



-V6s 2 (r-l)af, 



where 



T = V 



dV 



d 2 V 



The system (|12l) - (|T4| is non-autonomous since 

A 

Uch — 



(13) 
(14) 



(15) 



(16) 



and it is in general not closed since T does not depends 
a priori on the state variables x, y, s. Before to perform a 
detailed analysis of the stability of this dynamical system 
we rewrite it in terms of the observable quantities fi^ and 
7, and the new equations are given by 



Y = -37(2-7) + *(2-7)V37"* (17) 
= 3(1- 7)^(1 + 30^(1 -f^Kh, (18) 

s' = -V3s 2 (T- 1)V7^- (19) 

Equation (|17l) together with eq. (TT9")) encode the exact 
description of the dynamic evolution of the scalar field. 
In any case to find an exact solution it is a difficult task, 
and in order to proceed we consider two assumptions. 
First, we consider that the barotropic parameter of the 
scalar fluid is near to —1, therefore 7 1. Second, we 
consider a near flat potential that is s is approximately 
constant, say s w sq. 

In the limit 7 C 1 we obtain from (|18l) an approxi- 
mated equation with solution 



3(a+l) _|_ 



(20) 



fi^a) = 

satisfying fl^^a = 1) = fi^o, where 

13 = (l-a^Xx + ir^a, 



and x = A/ B. Thus, in the limit x ~ * is recovered the 
solution described by the expression (25) in Ref. [i~2j : 



fi„ 



1 



fi 



00 



1 a 



corresponding to stan- 



dard quintessence. From the solution (1201) and from equa- 
tion ([8]) we obtain the key formula 



Q + l 



(21) 



On the other hand, let us introduce the auxiliary func- 
tion fi = y/7. Assuming that is a monotonic function 



of the scale factor (in order to avoid that dfi^ /da = at 
any value a = a ), we obtain from (JTTJ) and (|18[) 



dfi 
dflj, 



-2) (3/i-V3«o. 



1) 



(22) 



6(Sl0 - 1)0^ (-(J- 2 +Lo c h 

where we have used the approximation s ~ sq — const.. 
Using the hypothesis 7 <§; 1 one is able to use the Taylor- 
expand the above differential equation around /1 = up 
to second order to obtain the approximated equation 



so 

dn^ Vs(uj ch + 

+ a 

{ijJch + l)(fi# - l)£l<f> 



■o(m) s 



(23) 



Substituting into the equation (|23[) the expression = 
w c /,(Q^) given by (l2~Tj) and integrating the resulting equa- 
tion with the initial condition /x = at fi^ = (which 
is true for the models we are considering here) we obtain 
the exact solution for a = 1 given by 



~7T 



(24) 



where 



yj fi {p^x - 1) + 1^/1 - fi* {Py/x + 1) 
^ (/Vx-i) V 3 /V* + 3 



(25) 

F(i7|m) and _E(fi|m) are the elliptic integral of the first 
and second kind respectively, with 



fi 



/fi, 



and 



(26) 



(27) 



/Vx + i 

Finally we obtain the expression 

1 + w = /i (fi^,) 2 . (28) 

Expression (|28| allows to obtain an exact solution for the 
dynamical evolutions of the barotropic index and using 
(|20|) we obtain to = w(a). We would like to note that our 



solutions have the right dust limit when A — ¥ [12| , this 
behavior is shown in Figures [1] (a) and (b). 

In fact, in the limit A — > (i.e., x ~^ 0), the deviation 
between our solution (|2"5|) and the solution (23) in [12J is 
given by the term 



s 2 (Cl 4 



I] 



where 



sn 2 



= [ — E [ sin 



(29) 



1) +2(fi* - l)tanh _1 



/fill + 



(fi - ^^(sin- 1 



(30) 
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O* a 
(a) (b) 

FIG. 1: (a) The value o/w/sq fs. assuming a nearly flat potential and u ~ — 1. 27ie dotted (red) line corresponds to the 
approximated solution discovered in flB] (equation (23) in flB]) and the dash-dotted (dark) line corresponds to the approximated 
solution, l|28[). for x = 0.1, presented here. 

(b) EoS parameter of the scalar field for the exponential potential with slope so = 0.8029. We have considered the initial 
conditions 7(1) = 0.1, 0^(1) = 0.7, s(l) — 0.8029. The scale factor is normalized to 1 at present. The continuous (green) line 
corresponds to the exact value of u>(a) for model with solely a scalar field; the dashed (blue) line corresponds to the exact value 
of w{a) for model that includes the Chaplygin gas; the dotted (red) line corresponds to the approximated solution found in flij/ 
(equation (23) in fl^l ) and the dash-dotted (dark) line corresponds to the approximated solution, (|28[) , presented here (assuming 
X = 0.1). 



But 



II. EXPONENTIAL POTENTIAL 



E($\m) 



1 — m sin 



2 a] 2 



dd 



[l - t 2 ] 2 [l - mi 2 ] 2 d*. 



(31) 



Thus, 




1 - sin 2 9] 2 d6 
(32) 



This means that h{Q,^) = 0. That is, in the limit A — » 0, 
our solution (|28l) and the solution (23) in 12] coincides. 

Particularly, in the Figure Q] (a) depict the value of 
w/sq vs. ri^ assuming a nearly flat potential and u ~ — 1. 
The dotted (red) line corresponds to the approximated 
solution discovered in [l2| (equation (23) in [12j]) and the 
dash-dotted (dark) line corresponds to the approximated 
solution, (|28j) . for x = 0.1, presented here. In Figure Q] 
(b) are displayed the EoS parameter of the scalar field 
for the exponential potential with constant slope. The 
scale factor is normalized to 1 at present. The continu- 
ous (green) line corresponds to the exact value of w(a) 
for model with solely a scalar field; the dashed (blue) line 
corresponds to the exact value of w(a) for model with 
the addition of the Chaplygin gas; the dotted (red) line 
corresponds to the approximated solution discovered in 
[l2| (equation (23) in [l2j|) and the dash-dotted (dark) 
line corresponds to the approximated solution, (|28[) , pre- 
sented here (assuming \ = 0.1). 



In order to motivate the analysis for a general (arbi- 
trary) potential, let us consider the simpler case of the 
exponential potential V{<j>) = Voe"^ [llHii]. 

In order to do the analysis from the dynamical systems 
viewpoint of the mixture of a scalar field with exponential 
potential and a Chaplygin gas, we need to consider the 
variables x,y,s defined in the previous section plus the 
new variable 



A 



Then, from the equations of motions (HJ),© and ^ we 
obtain the following autonomous system 



-3a; 



3 2 Xy 2 + [1 + x 2 - y 2 ] 



3az 2 



z' = 3z [l + a + x 2 - y 2 ] - 3z 2 - — 



x 2 — y 2 



(33) 



We note, that in the dust limit (A — > 0) the variable z 
becomes automatically zero, the last equation in the sys- 
tem (|33[) is satisfied identically (z = defines an invariant 
set) and the remaining equations (|33p corresponds to the 
usual exponential quintessence scenario [33j |. 

Now, it is convenient to express the observable mag- 
nitudes in terms of the phase space variables. These ob- 
servable magnitudes are the dimensionless Dark Energy 
density, fi^, given by ([TO]) ; the equation of state (EoS) 
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parameter of the dark energy given by 



LO = — 

P<i> 

the EoS of the Chaplygin gas 

Pch 
Pch 



1 1 

X — y z 
x 2 + y 2 ' 



LO c h 



1 



the total (effective) EoS given by 

Plot 2 2 

oJtot = =x - y 

Ptot 

and the deceleration parameter 



q 



(ay 



= -1 



[1 + x 2 



-y -z 



(34) 



(35) 



(36) 



(37) 



These expressions are valid not only at the fixed points 
but also they are valid in the whole phase space. Thus, 
we evaluate them at the fixed points in order to determine 
the type of solution that they represent. 

In table (UJ) we show the existence conditions for the 
real and physically meaningful (curves of) critical points 
of the autonomous system fl3"3")) associated to the expo- 
nential potential and also the values of the dark-energy 
density parameter £Ide, of the dark-energy EoS param- 
eter lode, of the total EoS parameter Lo to t and of the 
deceleration parameter q evaluated at them. 

Now, let us discuss in more details the stability con- 
ditions for the corresponding critical points. The critical 
point A associated to a matter dominated universe is a 
saddle point. Observe that the singular points B,C,D 
and F belong to the singular surface x 2 + y 2 = 1. In 
this case, both, the denominator and the numerator of 
Eq. (|33p vanish simultaneously. In this case the addi- 
tional eigenvalue due to the extra z-coordinate could be 
finite positive or infinite with undefined sign depending 
of how the point is approached. Thus, linear approxima- 
tion fails and we need to resort to numerical investiga- 
tion. The critical points B and C, corresponding to stiff 
solutions, are always unstable. B (C resp.) is a local 
source for A > -VE (A < v/6, resp.), otherwise they are 
saddles. This argument is supported by numerical stud- 
ies as shown in Figures [5] (a)-(d), for the values of the 
parameters in the typical intervals (that are determined 
by the bifurcation values). Critical points D and E are 
the usual quintessence solutions widely investigated in 
the literature (see for instance [Hj]). Then, the main dif- 
ference with respect to the results found in Ref. |33[ is 
that, for a > none of them can be a late time attrac- 
tor. For D this argument is based on numerical analysis, 
since D belongs to the singular surface x 2 + y 2 and in this 
case, both the numerator and the denominator of the the 
system (|33p vanish, and then, the linear approximation 
is not valid. In the case of the point E we support this 
result in the fact that there exist one eigenvalues with 
positive real part. This means that if we include a GCG 



in the background, we cannot get a stable solution dom- 
inated by the scalar field (D) or an scaling solution (E). 
If we restrict our attention to the invariant set z = 0, 
then we have that D is a stable one for A 2 < 3 and thus 
it can be the late time state of the universe. In this case 
the equation for z is vanished identically, and we do not 
require to include this variable in the analysis. D cor- 
responds to a dark-energy dominated universe, with a 
dark-energy EoS in the quintessence regime, which can 
be accelerating or not according to the A-value. Addi- 
tionally, this solution is free of instabilities. This point is 
quite important, since it is stable and possesses lode and 
q compatible with observations [33[. Point E is stable in 
the invariant set z = 0. It can attract the universe at late 
times (in case of a GCG behaving as dust), and it is free 
of instabilities. It has the advantage that the dark-energy 
density parameter lies in the interval < £Ide < 1, that 
is it can alleviate the coincidence problem, but it has 
the disadvantage that it is not accelerating and possesses 
lode = 0, which are not favored by observations [33| . 
However, let us remark that they are saddles for the full 
vector field. The solution F exists for A = 0. It represents 
a de Sitter solution which is stable but not asymptoti- 
cally stable (see the Appendix I A The curve of critical 
points G (which exists only for A = 0) is stable but not 
asymptotically stable, whereas, K is asymptotically sta- 
ble (see details of the center manifold calculations for 
both G and K in the Appendix [B]) . To finish this sec- 
tion let us proceed to the discussion of some numerical 
elaborations: 

• Fig. [2] (a) shows several orbits for the values of 
the parameters a — 0.5, A = 0. B and C are lo- 
cal sources; A is a saddle. F coincides with D 
(contained in the curve G) and it is stable but not 
asymptotically stable for /(0) > 0, otherwise it is a 
saddle (see Appendix I A 2|) . E does not exist. Any 
arbitrary point in the curve G (that exist only for 
A = 0) is stable, attracting an open set of orbits. 
The center manifold of K is stable. 

• In Fig. [2] (b) we presented several orbits in the 
phase space for a = 0.5, A = 1. The kinetic- 
dominated solution B and C are local sources; the 
matter dominated solution A and the scalar-field 
dominated solution D are saddles and the Chaply- 
gin gas dominated solutions (which also mimics a 
de Sitter solution) K is the attractor. D is a local 
attractor in the invariant set z — 0. 

• Fig. [5] (c) shows several orbits in the phase space 
for a — 0.5, A = 2. B and C are local sources; A, 
D and E are saddles and K is the attractor. E is 
a local attractor in the invariant set z = 0. 

• Finally, in Fig. [2] (d) we display several orbits in 
the phase space for a — 0.5, A = 3. B (the kinetic- 
dominated solution with Xx > 0) is a saddle point; 
C (the kinetic-dominated solution with Ax > 0) is 
the local source (unstable node); D does not exists; 



Cr. P. 
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yc 




Existence 
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LOtot 


9 


A 











always 





arbitrary 








i 
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B 


1 








always 


i 


1 


arbitrary 


1 


2 


c 


-1 








always 


i 


1 


arbitrary 


1 


2 




A 


V 1 _ "6" 


A 

u 


\ 2 c 


i 
i 




arbitrary 


1 4- x2 




E 


hi 

V 2 A 


/ 3 
V 2(A) a 





A 2 > 3 


A? 











2 


F 





1 





A = 


1 


-1 


arbitrary 


-1 


-1 


G 





yc 


1-2/c 2 


yl <l,y c /0,A = 


?/ C 2 


-1 


-1 


-1 


-1 


K 








1 


always 





arbitrary 


-1 


-1 


-1 



TABLE I: The real and physically meaningful critical points of the autonomous system S33\l associated to the exponential 
potential. Stability conditions and the values of the dark-energy density parameter Q.de, of the dark-energy EoS parameter 
lode, of the total EoS parameter uitot and of the deceleration parameter q. 



Cr. P. 


Eigenvalues 


Stability 


A 


-1, |, 3(1 + a) 


saddle 


B 


3, 3 — y/|A, undef. 


unstable 


C 


3, 3 + y/|A, undef. 


unstable 


D 


-3 + A 2 ,-i(6- A 2 ), undef. 


saddle 


E 


3(l + a),/T,/3+ 


saddle 


F 


—3, —3, undef. 


stable (see Appendix 1 A 11) 


G 


0,-3(1 + a), -3 


NH 


K 


-3,0, -3(1 + a) 


stable 



TABLE II: The real and physically meaningful critical points of the autonomous system f33]) associated to the exponential 
potential. Stability conditions, NH stands for non-hyperbolic. Observe that the critical points B, C, D and F belong to the 
singular surface x 2 + y 2 — 1. In this case both denominator and numerator of k33\) vanish simultaneously. In this case the 
additional eigenvalue due to the extra z-coordinate could be finite positive or infinite with undefined sign depending of how the 
point is approached. Thus, linear approximation fails and we need to resort to numerical investigation. 



A and E are saddles and K is the attractor. E is 
a local (spiral) attractor in the invariant set z = 0. 



PHASE-SPACE ANALYSIS WITHOUT 
POTENTIAL SPECIFICATION 



III. 



In order to transform the system (|12j) - (|T4|) to an au- 
tonomous one, first, it is necessary to determine a spe- 
cific potential form V{<j>) of the scalar field </>. However, 
using the above example as a motivation, one could al- 
ternatively handle the potential differentiations using the 
auxiliary variable s given by 



while keeping the potential still arbitrary . The next 
step is to introduce the function 



/ = s 2 (r - l) 



v"{<f>) V{4>f 



(39) 



v(cb) v(4>) 2 ' 

to be an arbitrary function of s. In fact, if / can be 
expressed as an explicit one- valued function of s, that is 
/ = f(s), then, it is possible to write a closed dynamical 
system for s and a set of normalized- variables. On the 
other hand, by giving f(s), non identically equal to zero, 
we obtain the expressions 



1 



dK, 



V(s) = e f: oTm dK V , 



(40) 
(41) 



(oo\ 1 The variable s is just a constant (s = A) for the exponential 
V((f>) ' ' potential V = Vbe" A f 
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a=0.5, h=0 ct=0.5, X=1 




(c) (d) 
FIG. 2: The phase space of the system 133\) for: 

(a) a — 0.5, A = 0. B and C are local sources; A is a saddle. F coincides with D (contained in the curve G) and it is stable but 
not asymptotically stable for /(0) > 0, otherwise it is a saddle (see Avvendix lA 2\) . E does not exist. Any arbitrary point in the 
curve G ( that exist only for X — 0) is stable, attracting an open set of orbits. The center manifold of K is stable, also K is. 

(b) a — 0.5, A — 1. B and C are local sources; A and D are saddles and K is the attractor. D is a local attractor in the 
invariant set z — 0. The solutions at the x-y plane correspond to those at Figure 2 in Js^j. 

(c) a = 0.5, A = 2. B and C are local sources; A, D and E are saddles and K is the attractor. E is a local attractor in the 
invariant set z — 0. The solutions at the x-y plane correspond to those at Figure 3 in J3^j. 

(d) a — 0.5, A = 3. B (the kinetic- dominated solution with Xx > 0) is a saddle point; C (the kinetic- dominated solution with 
Xx < 0) is the local source (unstable node); D does not exists; A and E are saddles and K is the attractor. E is a local (spiral) 
attractor in the invariant set z = 0. The solutions at the x-y plane corresponds to those at Figure 4 in JhUl] . 



where the integration constants satisfy V(so) = Vq, 
4>{sq) = 4>q 2 . Thus, it is possible to reconstruct the 
potential V by the elimination of s between (f40|) and 
(j¥T|) . For the usual cosmological cases the potential can 
be written explicitly, that is V — V((j>). The details of the 
method, coined "Method of /-devisers", were presented 
in [53]. This method has the significant advantage, that 
one can first perform the analysis for arbitrary potentials 



2 We would like to note that the requirement that / must be dif- 
ferent to zero exclude of this analysis the case of the exponential 
potential and for this reasons we studied the exponential poten- 
tial in the previous section separately. 



and then just substitute the desired forms, instead of re- 
peating the whole procedure for every distinct potential 
(see [53] and references therein). 

Then, from the equations of motions ([T]), ([2]) and (|3]) 
we result in the following autonomous system 



-3a; 



3 2 
2 SV 



3 



3z [l + a + 
— v6/(s)i. 



2 21 o 2 

x — y I — 3z — 



— 7/2:, 
2 

3qz 2 



x 2 — y 2 ' 



(42) 
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Cr. P. 


Xc 




Zc 


S c 


Existence 






Well 


Utot 


i 


A 











S c 


always 





arbitrary 








i 

2 


B(s*) 


1 








s* 


always 


1 


1 


arbitrary 


1 


2 


C(s*) 


-1 








s* 


always 


1 


1 


arbitrary 


1 


2 


Dl o*\ 
L>(S ) 


s* 

-I s - 




u 


s 


[s ) So 


i 
i 


_i + 3" 


arbitrary 




+ ~ 2~ 


E(s*) 


/ 3 1 


/ 3 





s* 


(s*) 2 > 3 















V 2 — 


V 2 (=*)" 






2 


F 





1 








always 


1 


-1 


arbitrary 


-1 


-1 


G 





Vc 


1-2/c 





y 2 < l,y c =£0 




-1 


-1 


-1 


-1 


K 








1 


S c 


always 





arbitrary 


-1 


-1 


-1 



TABLE III: The real and physically meaningful (curves of) critical points of the autonomous system Existence conditions 

and the values of the dark-energy density parameter Qde, of the dark-energy EoS parameter lode, of the total EoS parameter 
uitot and of the deceleration parameter q. We use the notation s* for the values of s = s* such that f(s*) — 0, and s c for 
denoting arbitrary values of s at equilibrium. 



Cr. P. 


Eigenvalues 


Stability 


A 


-f, |, 0, 3(1 + a) 


saddle 


B(s*) 


3, 3- y§s*, —s/Ef(s*), undef. 


unstable 


C(s*) 


3, 3+ y§s*, n/6/'(s*), undef. 


unstable 


D(s*) 


-3 + (s*) 2 , -|(6 - (s*) 2 ), -s*/'(0> undef - 


saddle 


E(s") 


3(1 + a), p- (s*), P + 


saddle 


F 


undef. , undef. , S + , S~ 


stable (see Appendix 1 A 21) 


G 


0,-3(l + a),-3,A+,A- 


NH, stable for /(0) > 0, y c > 0, saddle otherwise 


K 


-3, 0,0, -3(1 + a) 


NH (unstable) 



TABLE IV: The real and physically meaningful critical points of the autonomous system Stability conditions, NH stands 

for non-hyperbolic. We introduce the notations /3 ± (s*) = | ^—1 ± ^ 24( - £ ^^^y 7 ^ ^ J , S ± = — | ^1 ± \J~^~ |/(0) j , and A* = 

— | ^1 ± \Jl — |y 2 /(0)^ . Observe that the critical points B{s*),C{s*),D(s*) and F belong to the singular surface x 2 + y 2 — 1. 

In this case both denominator and numerator of \4<ty are vanished simultaneously. In this case the additional eigenvalue due 
to the extra z-coordinate could be finite positive or infinite with undefined sign depending of how the point is approached. For 
F there are two eigenvalues whose nature depends on the way that F is approached. For that reason they are undefined. Thus, 
linear approximation fails and we need to resort to numerical works. 



In table Hill we present the existence conditions for the 
real and physically meaningful (curves of) critical points 
of the autonomous system (|42[) . We use the notation s* 
for the values of s — s* such that f(s*) — 0, and s c for 
denoting arbitrary values of s at equilibrium. We display 
also the corresponding values of the dark-energy density 
parameter fi^, of the dark-energy EoS parameter of 
the EoS of Chaplygin gas ui c h, of the total EoS w tot and of 
the deceleration parameter q. In table IIVI are presented 
the stability conditions for the the critical points. 

Now, let us comment briefly on the stability and phys- 
ical interpretation of the critical points of . 

The curve of critical point A is always a saddle. It rep- 
resents cosmological solutions dominated by the Chaply- 
gin gas mimicking dust, this solution correlates with the 



transient matter dominated epoch of the universe. Ob- 
serve that the critical points B(s*), C(s*), D(s*) with s* 
such that f(s*) = and F belong to the singular surface 
x 2 + y 2 = 1. In this case both denominator and numera- 
tor of (j4"2|) are vanished simultaneously. In this case the 
additional eigenvalue due to the extra z-coordinate could 
be finite positive or infinite with undefined sign depend- 
ing of how the point is approached. For F there are two 
eigenvalues whose nature depends on the way that F is 
approached, for that reason they are undefined. 

For s* , the solutions B(s*) and C(s*) are past at- 
tractors or saddle points under the same conditions of 
the standard quintessence scenario [33] with the iden- 
tification s* = A (see table IIV[) . They represent so- 
lutions dominated by the kinetic energy of the scalar 
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field mimicking a stiff fluid. The solutions D(s*),E(s*) 
and F represents the scalar field dominated solution, 
the matter-scalar scaling solution and de Sitter solutions 
dominated by the potential energy of the scalar field, re- 
spectively. The main difference here with respect the 
standard quintessence scenario [33| is that D(s*),E(s*) 
are saddle points (we are considering a > 0). So, the 
standard quintessence solutions are not late time solu- 
tions in this scenario. For analyzing the important criti- 
cal point F, the linear approximation fails and we need to 
resort to numerical studies or include higher order terms 
in the analysis. In fact, following our approach in the 
Appendix IA 21 we find that actually F is the late-time 
attractor for /(0) > and a saddle for /(0) < 0. 

Combining expressions (1161) and (1351) we find that as 



t — > +oo, z — > 1 — x 1 — y 2 . Thus at late times we can 



approximate the system (|42[) by 



-:-!.-- + 3x 3 + \/^sy 2 



y = 



sxy + 3yx , 



-VEf(s)x, 



and the decoupled equation 



z' = 6zx 2 . 



(43) 



(44) 



If x — > x c as t —> +oo, then from equation (|4"4"]) 
follows that z increases without bound in contradiction 
with the boundedness of 1 — x 2 — y 2 . Thus, as time goes 
forward, x — > 0. Hence z — > 1 — y 2 where < y c < 
1. By calculating the critical points of the system (j43|) 
and analyzing their linear stability we find that the only 
candidates to be the late-time attractors are: 

• the curve G which have the following system of 
eigenvalues and eigenvectors: 



analysis) , this curve is actually non- hyperbolic and 
it is not normally hyperbolic anymore, thus we can- 
not obtain information about its stability looking 
at the linearization. This one is the main difference 
that appears when considering the extra direction 
z. 

• The other candidate is the curve K which have the 
following system of eigenvalues and eigenvectors 

-3, 0, 

7(3j»0,l|, {0,0,1}, {0,1,0} 

which is also normally hyperbolic (the center sub- 
space is the plane y-s is tangent to the line s — s c ). 

• Finally, both numerical simulations and analytical 
methods suggest that F (contained in the curve G) 
is an attractor for /(0) > 0, and for /(0) < 0, it is 
a saddle (see Appendix IA 2\i . 

The above heuristic reasoning suggest that the future 
attractor of the system (|4"2"j) is located at the curve G, 
which contains the especial point F, or it is located at 
the curve K . 

With the exception of the point F, which is dominated 
by a constant potential, G represents a class of solutions 
where neither the potential energy of the scalar field nor 
the Chaplygin gas dominates. On the other hand the 
curve K corresponds to purely Chaplygin gas dominated 
solutions (which also mimics a de Sitter solution). 

Indeed, using the Center Manifold Theory it can be 
proved that if the condition /(0) > 0, is satisfied, the 
curve of fixed points G is stable but not asymptotically 
stable. Applying the same procedure to the curve K we 
find that also the curve K is stable but not asymptoti- 
cally stable. The details of the calculation are presented 
in the Appendix [Cj 



0, 



{0,1,0}, {- 




where A ± = -f (l ± yj\ - |^/(0)J . Since the 

center subspace is tangent to the ?/-axis, follows 
that the curve is normally hyperbolic 3 . Then fol- 
lows the stability of G on the space (x, y, s). This ar- 
gument is not complete, since we have forget about 
what happens in the z-direction. In fact, in the 
general case (when the z-direction is included in the 



3 Recall that a set of non-isolated critical points is said to be nor- 
mally hyperbolic if the only eigenvalues with zero real parts are 
those whose corresponding eigenvectors are tangent to the set. 
In this case the stability of the set can be deduced by examining 
the signs of the remaining non- null eigenvalues (i.e., for a curve, 
in the remaining n — 1 directions) |67|| . 



A. An example: Cosh-like potential 

The cosh-like potential V{<j>) — Vq [cosh (£</>) — 1] has 
been widely studied in the literature (see for example 

m m, m n, m, us m ehzj ilo, ci m, iz3). F 0r 

this potential, 



/(«) = -£)(* + £)• 



(45) 



Observe that /(0) = \^ 2 > 0. This is the sufficient 
condition for the stability of the class of de Sitter 
solutions represented by the curve of critical points G. 
For this choice s* € {£, — £} ■ Al so f'( s ) = ~ s ) thus 
/'(£) = = f. For this choice the system ((42]) 

admits twelve (curves of) critical points denoted by 
A, B(-0, C(-0, D(0, D(-0,E(0, £(-£), 
F, G and K. 

To finish this section let us discuss some numerical 
simulations: 
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a=0.5, 5=1 a=0.5, 5=1 




FIG. 3: Projections of some orbits in the phase space of the system \4%ty f or a — 0.5, £ = 1. Observe that B(±£) and C(±£) are 
local sources and the point F located at the curve G is a local attractor. The curve K is stable, but not asymptotically stable. 
The scalar field- dominated solution _D(±£) are of saddle type, as well as the rest of the (curves of) fixed points. The scaling 
solutions _E(±£) do not exist. 



• In the Fig. [3] are displayed the projections of some 
orbits in the phase space of the system (|42)l for 
a = 0.5, £ = 1. -B(±£) and C(±£) are local sources 
and the point F located at the curve G is a local 
attractor. The curve K is stable, but not asymp- 
totically stable. The scalar field-dominated solu- 
tion D(±£) are of saddle type, as well as the rest of 
the (curves of) fixed points. The scaling solutions 
E(±£) do not exist. 

• In the Fig. [4] are displayed some projections of 
some orbits in the phase space of the system (|42p 
for a = 0.5, £ = 2. Note that B(±£) and C(±£) are 
local sources and the point F located at the curve G 
is a local attractor. The curve K is stable, but not 
asymptotically stable. The scalar field-dominated 
solution D(±£) and the scaling solution E(±£) are 
of saddle type, as well as the rest of the (curves of) 
fixed points. 

• In the Fig. [5] arc presented the projections of some 
orbits in the phase space of the system (|4"!2)) for 
a = 0.5, £ = 3. Observe that B(±f) and C(±£) are 
local sources and the point F located at the curve G 



is a local attractor. The curve K is stable, but not 
asymptotically stable. The scalar field-dominated 
solution D(±£) and the scaling solution E(±£) are 
of saddle type, as well as the rest of the (curves of) 
fixed points. 



Summarizing, for arbitrary potentials, using the nu- 
merical simulations to support our conjectures, and em- 
ploying analytical tools as the main proof, we have cor- 
roborated that F is a late-time attractor which is con- 
tained in the curve G. The curve of fixed points G and K 
are stable but not asymptotically stable. The numerical 
simulations suggest that the early time attractors are: 

• B(£) or C(-f) for < £ < V6; 

• £(-£) or C(0 for -V6 < £ < 0. 

The rest of the (curves of) critical points are saddle 
points. 
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0=0.5, 5=2 





a.0.5, 5=2 



a.0.5, 5=2 





FIG. 4: Projections of some orbits in the phase space of the system for a = 0.5, £ = 2. TVote that £?(±£) and C(±£) are 
local sources and the point F located at the curve G is a local attractor. The curve K is stable, but not asymptotically stable. 
The scalar field-dominated solution D(±£) and the scaling solution i?(±£) are of saddle type, as well as the rest of the (curves 
of) fixed points. 



IV. FINAL REMARKS 

In the present work we have studied the thawing dark 
energy scenarios with the Chaplygin gas as the other 
matter content of the universe and different kinds of self- 
interacting potentials for the scalar field. First, we ob- 
tain the exact solution for the cosmological equation of 
our model in terms of the elliptic function of the first 
and second kind and we obtain the right limit in the case 
of dust matter. Second, we used the cosmological model 
for the exponential potential for the scalar field and also 
we perform the dynamical systems analysis and we char- 
acterize the phase space of this system. We found the 
critical point of this system and also we studied the sta- 
bility of this point. The main characteristic of this phase 
space are: 

• The critical point A is always a saddle. It rep- 
resents cosmological solutions dominated by the 
Chaplygin gas mimicking dust, this solution cor- 
relates with the transient matter dominated epoch 
of the universe. 

• The critical points B and C corresponding to stiff 



solutions are always unstable. B (C resp.) is a local 
source for A > -\/6 (A < \/6, resp.), otherwise they 
are saddles. 

• The usual quintessence points D (scalar field- 
dominated solution) and E (the usual scaling scalar 
field matter solution) [33| cannot be the late time 
attractors due to the presence of a GCG with a > 
in the background. However, in the limit A — > 0, 
when the Chaplygin gas behaves as dust and the 
z- variable is not required, we recover the standard 
quintessence scenario [33j . This is a crucial differ- 
ence with respect to previous works in the litera- 
ture. 

• For a constant potential (exponential with A = 0), 
the solution F represent a dc Sitter solution which 
is stable but not asymptotically stable (see the Ap- 
pendix |XT]). 

• For a constant potential (exponential with A = 0), 
the curve of critical points G is stable but not 
asymptotically stable (see details in the Appendix 
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a=0.5, 5=3 a=0.5, 5=3 

.... -• .... B® 




FIG. 5: Projections of some orbits in the phase space of the system 14 2j) /or a = 0.5, £ = 3. B(±£) ana! C(±£) are toca/ sources 
and the point F located at the curve G is a local attractor. The curve K is stable, but not asymptotically stable. The scalar- 
field- dominated solution D(±£) and the scaling solution £J(±£) are of saddle type, as well as the rest of the (curves of) fixed 
points. 



• For an exponential potential (A ^ 0), if is asymp- 
totically stable (see details in the Appendix [Bjl . 

Finally for analyzing general potential V(4>) we used 
the "Method of /-devisers" and we obtained the critical 
point for this dynamical system and also we characterized 
the phase space and we studied the stability of the critical 
point. In this case, the principal characteristic of this 
phase space are: 

• For arbitrary potentials, the curve of critical points 
A, representing dust solutions is always a saddle. 
This solutions correlated with the transient matter 
dominated epoch of the universe. 

• Considering s* such that f(s*) = 0. For these s- 
values, the solutions B(s*) and C(s*) are past at- 
tractors or saddle points under the same conditions 
of the standard quintessence scenario [33| with the 
identification s* = A (see table ITVj) . 

• For the same s-values, the standard quintessence 
solutions D(s*),E(s*) are saddle points (we are 
considering a > 0). This is the main differ- 
ence here with respect the results for the standard 



quintessence scenario 33J. In the dust limit A — > 0, 
the standard quintessence scenario is recovered 33] 
since the z-variable is not required anymore in the 
dynamics. 

• For arbitrary potentials and provided /(0) > 0, F 
(contained in the curve G) is an attractor and for 
/(0) < 0, it is a saddle (see Appendix I A 2[) . 

• For arbitrary potentials and provided /(0) > 0, the 
curve of fixed points G is stable but not asymptot- 
ically stable (see appendix ICl) . 

• For arbitrary potentials, the curve K is stable but 
not asymptotically stable (see appendix [Cjl . 
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Appendix A: Stability analysis of the pure de Sitter 
solution 



2. Arbitrary potential 

For the stability analysis of the point F we introduce 
the local coordinates 



{x, l-y,z,s}= {x, y, z, s} e + 0(e) 2 



(A5) 



where e< 1, and y > 0, z > 0. Then the evolution of the 
linear perturbations is given by the equations 



./■ -- -3x+ \l -s, 



In this appendix we introduce local coordinates for an- 
alyzing the stability of the pure de Sitter solution given 
by the fixed point F. 



1. Exponential potential 

For the stability analysis of the point F we introduce 
the local coordinates 



{x, l-y,z} = {x, y,z}e + 0(e) 2 



(Al) 



where e<l, and y > 0, z > 0. 

Then the evolution of the linear perturbations is given 
by the equations 



y = -3y + -i, 



z = 3az 



3az 



(A2) 



The system (|A2I) admits the first integral zy a — c\ , where 
ci is an integration constant. Thus we can study the 
reduced system 



3 

V 1 = -3y+ -ciz/~ Q . 



(A3) 



The system (|A3j) admits the solution passing by (xo,yo) 
at time r = given by 



x(r) = x e 3r , 

y{T ) = ( Cl - e~ 3 ^ ( Cl - 2^ +1 ) ^ ~" 



(A4) 



where c\ = z y^, z = z(0). Observe that C\ -C 1 as 
far as yo and zq are small enough. We have (x,y, z) — » 
(0,2-3Ti (z y Q )^r ,2^fr (z y Q )^r) j as t V oo . Thus, 
for a given 5 > 0, and a > 0, it is possible to choose an 

initial state such that z yg<2x(|) 2 , which give a 
final state in a (^-neighborhood of the origin. This implies 
the stability, but not the asymptotic stability of F. 



V = -3y + -z, 



z = 3az — 



3az" 



2y 

-V6xf(0). 



(A6) 



The system (|A6[) admits the first integral zy a = c\ , where 
ci is an integration constant. Thus we can analysis the 
reduced system 



x = -:-!.r -.- \/ -s, 



9' = -3y + -ciy Q , 
s' = -V6xf(0). 



(A7) 



The system (IA7|) admits the solution passing by 
(xo,yo, So) at time r = given by 



x(t) 



^s e- 3T / 2 sinh(£r) 



- 3r / 2 i (2C cosher ) - 3sinh(Cr)) 

2l : 



(r) = 2-^r U-e-^+Vr ( Cl -2^ +1 ) 



a+l\ \ a + 1 



_ s e" 3r/2 (3sinh(CT) + 2f cosh(^T)) + 



+ 



2£ 

(4£ 2 - 9) e" 3r / 2 io sinh(C r) 



(A8) 



( 



where /3 = ±^9-12/(0). 

For the choice /3 2 < |, i.e., /(0) > 0, s) — » 

0,2-^tt (z ^)^tt ]2 ^tt (z y«)^r ;0 ) , as T ^ oo. 

For /3 = ±|, i.e., for /(0) = 0, (x,y,z,s) -t 
(^,2-^tt (z ^)^tt ,2^tt (z y^)^T ,s ) , as r -»■ oo. 

Combining the above arguments we obtain that for 
/(0) > 0, F is stable, but not asymptotically sta- 
ble. For (3 2 > §, i.e., /(0) ~< 0, the per- 
turbation values x and s diverges, and (y, z) — > 
(2-^tt (z y%)^ ,2^tt (z y£)^r) , as r -> +co. Thus, 
F is a saddle for /(0) < 0. 



14 



Appendix B: Center Manifold calculations for an 
scalar field with exponential potential 

For analyzing the stability of the curve of critical points 
G (which exists only for A = 0) we introduce the new 
coordinates 



y c (2y (y 2 c - a - l) + y c (-yl + z + 2a + l)) 



a + 1 



v 1 = x, 



V 2 = - 



(y 2 c - a - 1) (2yy c - y 2 c + z - l) 
a + 1 



(Bl) 



which are referred to an arbitrary point at G with coor- 
dinates (0,y c , 1 — Uc)- Applying the procedure, we find 
that the center manifold is given by the graph 



(y ' 1) 

(ui,vi,v 2 ) ■■ vi = 0(5),v 2 = — h — — h 



Ay 2 c (a + l) 



"l (y 2 c - a - 1) bu\ (yl-a-l) 
8y c 2 (a + l) 2 64y2( a + l)3 

+0(5),|i*i| <5}, 



(B2) 



where 6 is a small enough constant, and 0(5) denotes 
terms of five order in the vector norm. The dynamics on 
the center manifold is governed by the equation 

uj = 0(5). 

From this follows that G is stable but not asymptoti- 
cally stable. The center manifold of K is given by the 
approximated graph 



u 2 X 



(x,y,z) : x = —j= + 0(5), y = u, 
u 4 (a-l)X 2 2 



z = 1 - "x+y _ u2 + ° (5) ' ,ui < s ] ' (B3) 

where 0(5) denote terms of order 5 with respect the vec- 
tor norm. 

The dynamics on the center manifold of K is dictated 
by the gradient-like equation 



0(5). 



(B4) 



Since the origin is a degenerated minimum of the po- 
tential U(u) = follows the stability of K. 



Appendix C: Center manifold calculations for an 
scalar field with arbitrary potential 

For study the stability of G we resort the the Center 
Manifold Theory. Let us assume that < /(0) < |. 



Then, introducing the new variables 

y c (2y (y 2 c - a - l) + y c (-y 2 c + z + 2a + l)) 



m = 



Vl 



v 2 = 



a + 1 

(yl -a- I) (2yy c - y 2 c + z - l) 
a + 1 

6 - 8/(0)j/2 - V6) +4f(0)x 
2^6 - 8/(0)y c 2 ' 



V3 = 



sy/6-8f(0)y*-4f(0)x + VGs 
2^6 - 8/(0)j/2 



(CI) 



and applying the procedure, we find that the center man- 
ifold is given by the graph 

{(ui,Vi,V2,v 3 ) : vi = g(u\) + 0(5), 

« 2 = 0(5),«3 = 0(5),|«i|< < y} ) (C2) 

where 9 ( U1 ) = ^(^ + a Ul (a^y )(l> ;- a -i ) > § . g & 

small enough constant, and 0(5) denotes terms of five 
order in the vector norm. 

The dynamics on the center manifold is governed by 
the equation 

u[ = 0(5). 

Form this follows that G is stable but not asymptotically 
stable. 

For analyzing the case of complex eigenvalues (/(0) > 
|), we can introduce the new variables 



v 2 = v -^,v 3 



V 2 - V3 

2i ' 



for deriving the real Jordan form of the Jacobian. The 
procedure is straightforward and the result is the same. 

For analyzing the stability of the curve of critical points 
K we proceed as follows. 

Let us assume s c ^ 0. Introducing the new variables 

mi = s-s c -^xf(s c ), u 2 =y, Vl = ^xf(s c ), v 2 = z, 

(C3) 

and applying the procedure, we find that the center man- 
ifold is given by the graph 



(tti, tt 2 , vi,v 2 ) : v x = ^s c ul + ^uxulf(s c ) + 0(4), 



v 2 = -ul + 0(A),u\ + u% < 5} 



(C4) 



where 5 is a small enough constant, and 0(4) denotes 
terms of fourth order in the vector norm. But since v 2 = 
z > 0, it follows that u 2 should be zero. Thus, the center 
manifold of the origin is 

{(1*1,1*2,^1,1^2) : 1*2 = 0,ui = 0(4), 

v 2 = 0{±),u\ + u 2 2 <8} . (C5) 



15 



The dynamics on the center manifold is governed by From this fact follows the stability (but not the asymp- 
the equations totic stability) of the center manifold of the origin, thus, 

follow the the stability (but not asymptotic stability) of 
«i=0(4). (C6) K. 
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